Attractive instability of oppositely charged membranes induced by charge density 

fluctuations 
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We predict the conditions under which two oppositely charged membranes show a dynamic, at- 
tractive instability. Two layers with unequal charges of opposite sign can repel or be stable when in 
close proximity. However, dynamic charge density fluctuations can induce an attractive instability 
and thus facilitate fusion. We predict the dominant instability modes and timescales and show 
how these are controlled by the relative charge and membrane viscosities. These dynamic instabil- 
ities may be the precursors of membrane fusion in systems where artificial vesicles are engulfed by 
biological cells of opposite charge. 



The interactions between lipid bilayers control many 
biological processes such as membrane fusion, in which 
the domains separated by each bilayer merge In or- 
der to fuse, the bilayer must overcome an energetic bar- 
rier which can be reduced by an increase of the local 
bilayer curvature Q. Thus, small wave length undula- 
tions of bilayers at close proximity promote fusion. We 
show here that electrostatic interactions between oppo- 
sitely charged bilayers give rise to such fluctuations and 
predict the dynamics of instabilities related to the cou- 
pled distance and charge fluctuations. These dynamics 
may control the fusion of vesicles with oppositely charged 
cell membranes and our results indicate how the time 
scale for such fusion can be optimized as a function of 
the relative charges and viscosities. 

The interaction of similarly charged surfaces is repul- 
sive at the mean-field level and in the limit of relatively 
high salt concentration can be described by the Debye- 
Hiickel (DH) approximation to the Poisson-Boltzmann 
equation [3] ■ This repulsion between the surfaces inhibits 
surface undulations and thus opposes fusion, but fluctu- 
ations in the density of the charged lipids can sometimes 
lead to an instability, as was demonstrated theoretically 
and numerically y|. The case of oppositely charged 
bilayers, however, is qualitatively different because (as 
we show below) a system governed only by electrostatics 
is always unstable; the unstable mode may be related to 
charge fluctuations, height undulations, or a mixture of 
both, depending on the spacing between the layers, the 
charge densities and the viscosities. 

Recent experiments show that positively charged lipid- 
DNA complexes fuse with negatively charged cell mem- 
branes |(J. The subsequent release and transfection of 
DNA in such systems makes these complexes possible 
candidates for gene therapy. Understanding and control 
of the fusion process in this system of oppositely charged 
membranes is important so that one can minimize the 
time in which the complex is engulfed by the cell. In ad- 
dition, once the complex is inside the cell it is surrounded 
by a positively charged bilayer that originated in the cell 
membrane. Another fusion of these oppositely charged 
membranes is needed in order to release the DNA into 
the cytoplasm 0- Our predictions for the stability of 
the system as a function of the membrane charges and 



viscosities, may allow the optimization of the charge and 
viscosity of the lipid-DNA complex in order to enhance 
the fusion and transfection process. 

Parsegian and Gingell 8] solved the Poisson- 
Boltzmann equation using the DH approximation for 
two oppositely charged surfaces. We imagine that the 
charge density is controlled by mixing charged and neu- 
tral lipids; this is the case in both cellular membranes 
and in the lipid-DNA complexes. When the charge den- 
sities on the surfaces are different, there are counterions 
that must remain in the region between the surfaces in 
order to balance the electrostatic interaction. At close 
proximity, the pressure of these counterions leads to a 
net repulsion between the surfaces. For flexible surfaces, 
this short range repulsion inhibits undulations when the 
surfaces are at a small distance. Such undulations are 
needed in order for fusion to occur. Thus, within this 
simple picture in which the surface charges are uniform, 
oppositely charged surfaces do not necessarily attract and 
will not necessarily fuse. In this paper, we consider the 
effects of lateral fluctuations of the charges in each mem- 
brane and show that at relatively short distances, these 
fluctuations give rise to a dynamic instability that pro- 
motes fusion. The instability involves a coupled mode of 
the local charge densities and the distance between the 
bilayers. 

In the DH limit, applicable to physiological salt con- 
centrations, the screening length that governs the inter- 
action between the bilayers is where k 2 = 8 ^y 3 , 
n s is the salt concentration, Q s is the charge per lipid and 
e is the dielectric constant of water. The charge density 
of each membrane i = 1,2 is cr^e/a, where a is the area 
of a lipid and <7j is a dimensionless charge density. The 
electrostatic free energy of the system per unit area is 
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where h is the dimensionless distance between the mem- 
branes in units of The last two terms in Eq. ^ 
account for the electrostatic repulsion between lipids in 
the same bilayer. 
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FIG. 1: The projections of the unstable eigenmodes of the 
system on the modes predicted for the limiting cases, (a) We 
use / = / e i and <7i = 0.8 to calculate the energy stability 
matrix, S, at each point on the h x 02 surface. Let u be the 
eigenvector with the negative eigenvalue of S. In each region 
the color code represents the projection u • Vi, with v% is vu, 
v ai , 1V2 or (defined in the text), according to the label 
on that region. Below the black-grey dashed line fhh > 
so a a system with fixed constant charge density profile is 
stable there, (b) The same, but for / = / e i + /ent- Below the 
dashed gray line all the eigenvalues are positive, and there is 
no unstable mode. We choose u there as the eigenvector with 
the smallest eigenvalue, (c) The same as in (a), but for the 
unstable eigenvector of the dynamical matrix, D, (see Eq. |^. 
(d) The same as in (b), but for the eigenvectors of D. 

We assume that the average distance between the bi- 
layers is imposed by external constraints, such as lateral 
or transversal pressure, binding of trans-membranal pro- 
teins, or the system geometry, as in the case where the 
cell membrane encloses the lipid/DN A complex 7J. Sim- 
ilarly, the overall charge on each membrane is conserved. 
We therefore focus only on local changes in the inter- 
membrane spacing and charge density. 

As a first step, we consider only the electrostatic free 
energy, that is, we let the total free energy, / = / e j. In 
the calculations below, we assume o\ > ~a 2 > 0. The 
thermodynamic stability (see below for the dynamics) of 
the system depends on the eigenvalues of the Hessian 
matrix S of the second variations of the free energy /, 

(fhh fhai fh(T2 \ 
fh<j\ fa\<j\ fa\<J2 I ■ (2) 
fha 2 JO 102 f&2&2 ) 

The determinant \S\ = (fffy) s^Mh) 1S a l wa Y s negative 
for oppositely charged layers. In fact, S has one negative 
eigenvalue Ai, which corresponds to the unstable eigen- 
mode Vi of the system. We estimate this unstable mode 
analytically for three limiting cases: 

(1) For /i » 1 we have A;,, = fffrcrio^ exp(— h) and 



Vh = (1,0,0). This represents a mode where there are 
surface undulations but with a constant charge density. 
This constant density is the result of the repulsion be- 
tween similar charges on the same bilayer, which inhibits 
charge density fluctuations. In this case, the instability 
is due to the "naive" picture of attraction of oppositely 
(but uniformly) charged surfaces. 

(2) When h <C 1 and h <C o\ + 02 the non-diagonal 
terms of S vanish and h, o\ and <j 2 decouple. The un- 
stable mode (under our assumption of \(J%\ < is 
v a2 = (0,0,1) with X a2 = 4 ™ a2 gi °^ ff2 h (in the case of 
I £72 1 > I £7i I the unstable mode is v ai = (0,1,0)). The 
associated density fluctuations tend to lead to regions 
where locally o\ ~ —a 2 @- This is the situation in 
which the largest attractions occur since there is no need 
for countcrions to balance the electrostatic interaction. 
Thus, the counterions may be depleted from these regions 
and the associated entropic repulsion between the lay- 
ers vanishes. This situation, where the interaction takes 
the form of a nearly bare Coulomb attraction between 
surfaces, is energetically favorable. There are fluctua- 
tion only in £72 since the repulsive interactions between 
charges in the same bilayer inhibits charge fluctuations 
in the more highly charged bilayer (a 2 2 < ci 2 ) and the 
charge fluctuations are therefore limited to the bilayer 
with the lower charge density. 

(3) When <j\ + a 2 <C h <C 1 the charge densities 
are almost equal and opposite, o\ ~ —£72, and the in- 
teraction between the layers takes the form of an un- 
screened Coulomb attraction. Charge density varia- 
tions serve to increase the instability. The eigenmode 
v m = I (%/2, — 1, l) is an equal mixture of height undula- 
tions and in-phase charge density fluctuations. The en- 
ergy associated with this instability is given by the eigen- 

value A m = — v % \o\ — 0-2 1 • In order to estimate the 
contribution of the charge fluctuations to this instability, 
we compare its energy scale, given by A m , with the energy 
scale for the same instability of a system with uniform, 
local constant charge densities, given by fhh- To zeroth 
order in <ti + a 2 we have ^ — hm , h o^/^. ~ h <C 1. 

A m v2cosh J (/i/2) 

This means that charge fluctuations, when they are al- 
lowed, significantly increase the unstable eigenvalue and 
thus speed up the dynamical attractive instability. 

In order to more accurately determine the regions of 
validity of these limiting cases, we calculated numerically 
the unstable eigenmodes of the stability matrix, S 1 , and 
their projections on each of the "ideal" modes found an- 
alytically in the limiting cases discussed above. As seen 
in Fig. nji, for h < 1 there is wide ranges of values of 
£72 where the mixed mode v m is dominant. The charge 
density modes v ai and v a2 are dominant at h < 0.25 and 
|l + £7 2 /£7i| >0.3. 

In addition to the electrostatic energy, charge fluctua- 
tions caused by the demixing of the charged and neutral 
lipids also modify the local entropy of each bilayer, given 
by fi = [ct logc, + (1 - Ci ) log(l - a)] (for i = 1, 2), 
where Cj = UiejQi is the number density of the charged 
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lipids in bilayer i, and Qi is the charge of a single lipid. 
The entropy will tend to remix the two species and may 
inhibit to some degree or even eliminate the charge in- 
stabilities that are promoted by the electrostatic interac- 
tions. On the other hand, attractive interactions, such as 
Van der Waals attraction between similar lipids, promote 
the demixing charge fluctuations and enhance the insta- 
bility. Here we consider the "worst case" where there are 
no such attractions. 

We also include the Helfrich repulsion between the 
bilayers, induced by the entropy of the bilayer undula- 
tions 10]. Its contribution to the free energy is /h = 
cu(nkBT) 2 (ki 1 +k2 1 )h~ 2 , where ch — 0.116 is a univer- 
sal number, and fej is the bending modulus of bilayer i. 
The contribution of the entropy of the charged lipids and 
of the undulations to the free energy is / en t = fi + f 2+ fa- 
in Fig. we present results for the eigenmodes of the 
stability matrix, S, for the total free energy / — / c i+/ e nt- 
The smallest eigenvalue is positive in the regions where 
the v rTl and v a2 modes arc dominant, and the system is 
stable there. Nevertheless, there is a finite range of er 2 
values, for which the mixed mode v m is dominant, where 
the system is unstable even for values of h as small as 
0.4. Thus, the entropy modifies, but does not eliminate 
entirely the predicted instabilities. 

In order to test this prediction one can measure the 
stability of multilayer systems made of alternating layers 
of positive and negative charge. Systems that contain 
random mixtures of cationic and anionic lipids were al- 
ready examined [lll | . However, for the instabilities to be 
observed, the cationic and anionic lipids must phase sepa- 
rate (due to e.g., different chemistry of the polar heads or 
tails). The spacing between the layers can be determined 
by the fraction of water in the mixture. We predict that 
the system will be stable when the interlayer distance h is 
small, but it will become unstable when h is larger than 
a critical value h c as shown in Fig. ^ which depends on 
the charge densities of the bilayers. 

We now consider the dynamical response of the system 
in order to predict not just the regions of possible stabil- 
ity, but also the time scales associated with the distance 
and charge instabilities. These timescales and their de- 
pendence on the charge and the viscosities are important 
in controlling the fusion processes in the experimental 
systems described above. We therefore consider the dy- 
namical response of the system to a small perturbation 
with a long wavelength 2irq~ 1 ^> h, where h = n~ 1 h is 
the local instantaneous distance between the membranes, 
and q is the wave-number of the perturbation. 

The hydrodynamics of the water layer between the bi- 
layers in the approximation of long wavelength pertur- 
bations is given, to first order in qh, by the lubrication 
approximation 1 2] : 
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where r\ is the water viscosity and 7 is the membrane 
surface tension, which we can neglect if the wave number 
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FIG. 2: The time scale r for the growth of the instability for 
different values of C,\. t is given by inverse of the negative 
eigenvalue of the dynamical matrix, D, calculated for / = 
/ci + /ent, h = 1 and (2 = fO -7 erg s/cm 2 . We keep o\ — 0.8 
and vary 02 along the x-axis. The dotted line is the time scale 
for a system with a uniform charge density, r diverges when 
the system is stable. 



q is small enough (see discussion at the end of the letter). 

The dynamics of the charged lipids in the membrane 
is governed by the Smoluchowski equation 0] 



dt 



d_l_ 

dx Ci 



' dx 



df_ 

da. 



(4) 



where fii (i = 1,2) is the chemical potential of a charged 
lipids in bilayer i, Ci is the number density of the charged 
lipids and Q is the two-dimensional viscosity in the bi- 
layer. 

We consider a perturbation of the type h(x, t) = ho + 
5h(t) sm(q-x), <Ti(x, t) = <Ji + 5ai(t) sin(q-x). From Eqs. 
[3] and 0] we obtain the dynamics to first order in Sh and 
5(Ji\ 




(5) 



where the derivatives in S are evaluated at h = ho and 
From Eq. [5] we obtain characteristic time scales for 



e^a — 2 
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the dynamics h, <J\ and oi'- — 3r/K 
rtrr^Q 2 - For the numerical calculation of the time 
scales we use the values: o\ = 0.8, Qi — —Q2 = 2e, 
r] = O.Of erg s/cm 3 , £i = £2 = 10~ 7 erg s/cm 2 [T^ . 
T = 298 K, ki > k 2 = 20k B T, a = 10" 14 cm 2 and 
k = 10 7 cm -1 . With these values we have r CTi = 35t/j. 



cm 

If we consider e = 47re e W ater = 
and q = 10 5 cm -1 we have Th 



3.5 • 10 20 e 2 erg" 
= 1.7- 10" 6 s. 
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The time scale for the growth of the instability is the 
inverse of the negative eigenvalue of the dynamical ma- 
trix, D (in the case that one exists) . We have calculated 
the eigenmodes and time scales for the two cases we con- 
sidered: (i) electrostatics alone, / = f c \ (ii) electrostatics 
and entropy, / = / C i + / C nt- In Figs. ^ and^i we present 
the regions where the unstable mode corresponds to one 
of the modes we found above for the limiting cases. The 
figures demonstrate that our results, obtained through 
analytical approximation, are valid for large regions in 
the parameter space. 

When we take into account the entropy (Fig. (IJi) the 
mode v CT2 almost vanishes and the regions in which the 
modes v a2 and v ai were unstable in the absence of en- 
tropy, become stable. However, there is still a range of 
charge densities for which the system is destabilized by 
charge fluctuations. This range is represented in Fig. ^1 
by the area between the two stability lines. It is mostly 
dominated by the mixed mode v m that couples the charge 
density and distance fluctuations. This mixed mode re- 
sults in a time scale for the attractive instability that 
is significantly faster than the time scale predicted from 
the mean- field attraction (i.e., in the absence of charge 
density fluctuations). 

In Fig. [3 we present the dependence of the time scale 
t for the instability on the intra- membrane viscosity £i. 
As expected, we see that if £i > (2 the time scales in the 
region where |(72| > \<Ji\, in which v ai mode is dominant, 
are much slower than the time scales in the opposite re- 
gion, where v a2 mode is dominant, and vice versa. An 
increase in the viscosity £1 shifts the maximal instability 
(and shortest time scale) to lower values of the charge 
ratio. This prediction indicates how one might optimize 
of the charge and viscosity values in the fusion and trans- 



fection experiments. 

The electrostatic free energy / c i promotes that growth 
of perturbations with any wave number q. These insta- 
bilities in are inhibited by the surface tension 7,; and 
the bending modulus ki of each membrane labelled by 
i = 1,2. For simplicity, we assume here that = k <C k\ 
and 72 = 7 <C 71, so that variations of the inter- 
membrane distance, h, only arise from undulations of 
the more flexible bilayer, denoted as bilayer 2. In this 
case, the only contribution of these terms to the sta- 
bility matrix, S, (defined for the electrostatic-only case, 
/ = /el) will be in the entry S n = U\ h h + l<f + §<7 4 , and 

I o| _ 167T 2 e 4 / 47rg-KT 2 e 2 , 2 , k A 
PI — £ 2 K 2 a 4 \e K a 2 sinh(/i) + 11 + 2^ )' 

Let g max be the wave number at which the system is 
marginally stable: \S\ = 0. The system is unstable for all 
perturbations with longer wavelengths, q < <7 max - From 
Eq. [SI one sees that the time scale for the increase of a 
perturbation is proportional to q~ 2 . Thus the dominant 
instability of the system will be with a wave number q ~ 

?max- 

At a small inter-membrane separation h <C 1 we have 
<Zmax ~ h~* . As h — > our theory predicts the growth of 
unstable modes with large wave vectors. The related un- 
dulations will grow until the bilayers are at a microscopic 
distance, where the electrostatic instability is balanced 
by short-range hydration repulsion. Thus, there will be 
regions where the bilayers are at close proximity and have 
high local curvature. In these regions the barrier for fu- 
sion is reduced considerably |2j, and the probability for 
fusion is greatly increased. 
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